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A current threat to the marine ecosystem is the high level of solar ultraviolet radiation (UV). Large whales 
have recently been shown to suffer sun-induced skin damage from continuous UV exposure. Genotoxic 
consequences of such exposure remain unknown for these long-lived marine species, as does their capacity 
to counteract UV- induced insults. We show that UV exposure induces mitochondrial DNA damage in the 
skin of seasonally sympatric fin, sperm, and blue whales and that this damage accumulates with age. 
However, counteractive molecular mechanisms are markedly different between species. For example, sperm 
whales, a species that remains for long periods at the sea surface, activate genotoxic stress pathways in 
response to UV exposure whereas the paler blue whale relies on increased pigmentation as the season 
progresses. Our study also shows that whales can modulate their responses to fluctuating levels of UV, and 
that different evolutionary constraints may have shaped their response strategies. 

Although the Montreal Protocol, (which in 1987 banned the use of ozone depleting substances) has been 
central to decelerating the loss of ozone^ substances released during the nineties continue to destroy the 
ozone today^'^. The resulting large amount of solar ultraviolet (UV) radiation reaching the biosphere 
represents a significant threat to our ecosystems and there is now evidence that, as is the case in humans^ wild 
species such as amphibians, fish, and invertebrates suffer from sun-induced skin damage'*"^. Until recently, effects 
of UV on large marine mammals, which due to their life history and physiological constraints are unable to avoid 
continuous exposure to the sun, have been ignored^. Nevertheless, due to their long life expectancy and their 
extended distribution that spans all oceanic latitudes, cetaceans could reflect UV variation across large spatial and 
temporal scales and thus be considered as "UV-barometers of the ocean". Although it has been demonstrated 
recently that UV-induced acute damage is widespread and significant in whales^, nothing is known about the 
genetic nature of this damage, or the strategies used by these animals to resolve or counteract UV-induced lesions. 
Following initial observations which revealed that having fewer melanocytes led to more UV-related skin lesions 
and fewer apoptotic cells amongst whale species, it was suggested that darker pigmentation confers cellular 
protection from sun irradiation and plays a role in the elimination of potentially precancerous cells in whales^. 
Whether such a role has an effect at the molecular level remains unknown. Considering that UV-induced DNA 
mutations can contribute to the development of skin cancer in humans^, it is relevant to record the occurrence and 
magnitude of UV-induced DNA damage in whales. One particularly successful biomarker of UV-exposure in 
human skin is mitochondrial DNA (mtDNA)^°. Because mtDNA has a higher rate of mutation than nuclear DNA 
and a reduced capacity to repair damage^^, UV-induced mtDNA lesions accumulate throughout the life of an 
individual and thus offer a reliable biomarker for cumulative exposure to UV^°'^\ Mitochondria are the predom- 
inant site for the formation of cellular reactive oxygen species, which tend to be exacerbated following UV 
exposure in human skin causing mtDNA damage in the form of strand breaks and deletions ^\ Furthermore, 
accumulated mtDNA damage has been associated with human skin cancer^ \ We sought to test two hypotheses: 
first, that whales accumulate mtDNA damage in their skin as they age, and second, that pigmentation confers 
protection against UV-induced mtDNA damage. We screened mtDNA from skin of three whale species (fm, 
sperm, and blue whales) searching for UV-induced damage using quantitative long- PGR assays and investigated 
their association with minimum age and with "individual pigmentation indices" including melanocyte abund- 
ance and melanin density. We next sought to investigate the capacity of whales to respond to UV and study 
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interspecies differences in the mechanisms used in this process. For 
this, we tested three hypotheses: first, that UV-induced damage leads 
to changes in the expression of genes involved in genotoxic stress 
pathways; second, that species with different skin colour and sea 
surface behaviour (time spent at the sea surface) use different cellular 
mechanisms to protect themselves from UV exposure; and finally, 
that whales modulate UV-counteractive molecular mechanisms in 
response to the seasonal increase in UV levels, as is known to occur in 
humans^^. To test these hypotheses, we used quantitative PGR assays 
to monitor the expression of carefully selected genes HSP70, 
encoding heat shock protein 70, an indicator of cell stress^^'^^, tumour 
protein P53 (P53), a central transcriptional factor activated by stres- 
sors such as UV and involved in cell cycle arrest, DNA repair and 
apoptosis^^'^^ and KIN17 (KIN), a cell cycle control protein 
up-regulated by UV^^'^° and constructed linear mixed effect 
models To fully explore the role that cetacean skin pigmentation 
has in shaping protection against UV exposure, we also measured 
transcription levels of the tyrosinase gene (TYR), a key player of 
melanogenesis^^'^l 

Results 

Do whales accumulate mtDNA damage in their skin as they age, 
and does pigmentation confer protection against UV-induced 
mtDNA damage? Using quantitative real time PGR methodology 
(qPGR), we screened for mtDNA damage from the skin of three 
whale species (fin, sperm, and blue whales; Fig. lA) with different 
skin colour (Fig. lA-G) and sea surface behaviour. The level of blue 
whale mtDNA damage was inversely predicted by melanin density of 
individual samples (tn = —2.65, p = 0.02; see details in Table SIA; 
Fig. 2G), which suggests that greater pigmented individuals 
accumulate lower levels of damage. In addition, the abundance of 
melanin was also inversely related to the level of microscopic lesions 
(cytoplasmic vacuolation and intracellular oedema, methodology 
described previously^)) for whales across the three species studied 
(t95 = -2.11, p = 0.04 and tg^ = -4.19, p<0.0001 respectively; see 
details in Table S2; Fig. 2 A and B). Furthermore, as predicted under 
the assumption that mtDNA damage accumulates over time within 
the cell (aided by compromised mtDNA repair mechanisms ^°), the 
amount of whale skin mtDNA lesions increased with age for blue 
whales (tn = 2.24, p = 0.05; see details in Table SIA; Fig. 2E). 

Melanin abundance was not only dependent on the quantity of 
melanocytes (Fi,4o = 19.83, p<0.0001; Fig. IG) but also on the 



transcriptional activity of TYR and P53 (¥2,40 = 6.43, p = 0.004 
and Fi 40 = 4.06, p = 0.05, respectively), genes known to be involved 
in the process of melanogenesis^^"^^. The abundance of melanocytes^ 
and of melanin was highest for fin whales, the darkest of the species 
studied here (tgs = 3.99, p = 0.0001 and tgs = 3.42, p<0.001 when 
compared with sperm and blue whale, respectively; see details in 
Table S2; Fig. IB and Fig. SIA). Intriguingly, both blue and sperm 
whales shared comparable levels of melanin (tg^ = —0.93, p = 0.35; 
see details in Table S2), although melanocytes were more prevalent in 
sperm whales than in blue whales^, suggesting that melanin produc- 
tion capacity is restricted in sperm whales compared to blue whales. 

The capacity of whales to respond to UV and the interspecies dif- 
ferences in the molecular mechanisms used in this process. To 

accomplish this we tested three hypotheses. Firstly, we predicted 
that UV-induced damage would lead to changes in the expression 
of genes involved in genotoxic stress pathways in the studied whale 
species. Secondly, species with different skin colour and sea surface 
behaviour (time spent at the sea surface) might use different cellular 
mechanisms to protect themselves from UV exposure. Thirdly, 
whales would expectedly modulate UV-counteractive molecular 
mechanisms in response to the seasonal increase in UV levels. 

As part of the process in addressing the first two hypotheses, the 
level of HSP70, an indicator of cellular stress^^'^^, was measured using 
qPGR. Transcription of HSP70 was significantly higher in sperm 
whales when compared with fin and blue whales {t^o = 4.57, p = 
0.0001 and t^o = 3.75, p<0.001 respectively; see details in Table S3B; 
Fig. SIB). In addition, HSP70 expression and mtDNA lesions were 
found to be inversely related across species (ti5 = —4.26, p<0.001, 
see details in Table SIB; Fig. 2D). To complete the process of addres- 
sing the two hypotheses, the expression of KIN, a cell cycle control 
protein up-regulated by UV was determined^^'^°. The results show 
that KIAT expression was also higher in sperm whales {t^^ = 3.12, p = 
0.004 and tss = 3.95, p<0.0005; see details in Table S3A; Fig. SIB). 

To address the third hypothesis of seasonal changes in stress res- 
ponse, we found that when studying temporal changes in gene 
expression levels, the transcription levels of P53 and HSP70 formed 
an ascending curve between February and May, peaking in March/ 
April (LR = 11.93, p<0.01 and LR = 8.96, p = 0.03, respectively; see 
details in Table S3; Fig. 3A). Following the seasonal increase in UV 
radiation that reached the Gulf of Galifornia (Fig. 3B), blue whales 
increased melanocyte abundance (Fi^47 = 3.07, p = 0.04; see details 
in Table S4; Fig. 3G), suggesting that individuals of this species 
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Figure 1 | The contrasting skin colours of the three studied species. (A) Photograph showing, from top to bottom, a blue whale (pale grey skin colour, 
the lightest species), a sperm whale (dark grey skin colour) and a fin whale (black skin colour, the darkest species). This photograph was taken by DG. 
(B) Differences in density of melanocytes (thick bars) and melanin (thin bars) amongst the three studied species (n = 53, n = 17, n = 45 for blue, sperm 
and fin whales respectively). (C) Association between melanin abundance and melanocyte counts in whales. Grey dots correspond to blue whales, 
black dots to fin whales and crosses to sperm whales. The counting area was determined as previously described (see^). Briefly, the number of melanocytes 
per 100 arbitrary units was determined in triplicate along the entire epidermal ridge. The number of epidermal ridges to count was established a priori 
based on a cumulative curve*^. 
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Figure 2 | Melanin abundance, HSP70 expression and age influence sensitivity to UV-induced damages in whales. Relationship between whale melanin 
abundance and skin lesions: (A) intracellular oedema (n = 39 and n = 66 for absence and presence, respectively) (B) cytoplasmic vacuolation (n = 9, 
n = 23, n = 41, n = 35 for vacuolation levels of 0, 1, 2 and 3, respectively; see^), and (C) blue whale mitochondrial DNA lesions. (D) Inverse 
correlation between mitochondrial DNA (mtDNA) lesion density and HSP70 expression levels (ACt). (E) Direct correlation between mtDNA lesion 
density and individual minimum age of blue whales (calculated by taking into account the first year of observation reported for a particular individual in 
the Gulf of California; individuals with a minimum age of 1 were excluded; individuals that were observed in the Gulf of California on the year they were 
born, and as such their exact age was known, were included). Bars = ± 95% CI. 



are capable of modulating their degree of pigmentation. Although 
statistically non- significant, melanin density appeared to follow the 
same increasing trend over months suggesting tanning ability (see 
details in Table S4; Fig. 3C). Intriguingly, we did not observe ana- 
logous changes in the pigmentation of the (comparatively darker) fm 
whale (Fig. 3D). 

Another, non-exclusive, hypothesis for the observed differences in 
pigmentation plasticity amongst species might entail their distinct 
migratory behaviour. While fin whales are year-long residents of the 
Gulf of California^^, blue whales migrate annually from higher to 
lower latitudes^^, where levels and intensity of UV are greater^^. 
Consequently, when blue whales arrive at the Gulf of California they 
are exposed suddenly to higher levels of UV. It is possible that the 
higher prevalence of microscopic lesions and mtDNA lesions that 
occur at the beginning of the season (Fig. 4 A, Table SI A) reflects the 
time needed for UV acclimatization to occur^°, especially as blue 
whale pigmentation seems to then increase gradually (Fig. 4B). 

Discussion 

Those blue whales with higher levels of melanin showed lower levels 
of mtDNA damage; analogous correlations have been recorded 
recently in humans, strongly suggesting that melanin confers pro- 
tection against mtDNA UV-induced damage^ \ Those whales with 
higher melanin also showed a reduced level of microscopic lesions. 
Our results thus offer further evidence that melanisation protects 
animals such as whales against UV exposure, at both the molecular 



and cellular levels^. The level of mtDNA damage also increased with 
age in those whales for which information was available (the blue 
whales), which is also observed in humans and other species^°, and 
could result in decreased bioenergetic capacity with age. 

The transcription of HSP70 was significantly higher in sperm 
whales compared to blue and fin whales, which could suggest a 
potential inhibitory effect of HSP70 on melanin production, as has 
been observed in mice that over-express HSP70^^. Paradoxically, 
HSP70 appears to have a geno -protective effect for cells exposed to 
UV. Evidence for this assumption is supported by the observation 
that those whales with higher levels of HSP70 showed decreased 
mtDNA lesions across species. Expression of KIN was also higher 
in the sperm whales. Taken together, transcription patterns of HSP70 
and iCJAT might reflect the markedly dissimilar sea-surface behaviour 
of sperm whale compared to the other two species. For instance, 
foraging dive intervals at the surface are approximately five times 
longer for sperm whales than for blue and fin whales^^'^^, and during 
socialization, sperm whales can remain at the surface for up to six 
hours at a time^^, thus increasing exposure to damaging UV. In 
humans and laboratory animals, levels of expression of repair genes 
such as HSP70 and KIN increase in a time-dependent manner fol- 
lowing UV irradiation and up-regulation of HSP70 and iCJATcan 
be observed between 6 h and 8 h post UV irradiation^^. In this sense, 
the recorded expression levels suggest that sperm whales activate 
genotoxic stress pathways that involve over- expression of HSP70 
and KIN in response to long and persistent exposure to UV. 
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Figure 3 | Monthly fluctuations in gene expression and pigmentation levels follow seasonal variation in UV. (A) Monthly differences in mean 
expression levels of P53 and HSP70 genes in blue and fin whales (data pooled; expressed as ACt, y axis inverted; n = 11,13,12 and 6 for February, March, 
April and May, respectively). Sperm whales were not included as they were sampled exclusively in May. (B) UV index recorded between January and June 
over the Gulf of California, Mexico (data average records for 26°-28°N and 109°-112°W) for the years 2007 (red), 2008 (blue) and 2009 (green). 
Calculation (a simply function of total column ozone and the solar zenith angle) was conducted under local noon and clear sky conditions and does not 
consider cloud or aerosol effects. Observation years extend from 1979-2010 (32-year running average shown by a thick black line). The lower and upper 
thin black lines show the minimum and maximum value observed, respectively. The grey shading shows the probability distribution function (i.e., 80% of 
the observations are within the light grey shading, while 40% are within the dark shading). Plot obtained using total ozone observations from the Total 
Ozone Mapping Spectrometer (TOMS) and Ozone Monitoring Instrument (OMI). This figure was constructed by Eric Nash and Paul Newman from 
NASA Goddard Space Flight Center. (C) Monthly variation of blue whale melanocyte and melanin abundance during 2007 (n = 3, 13, 7 and 3 for 
February, March, April and May, respectively). (D) Monthly variation of fin whale melanocyte and melanin abundance during 2008 (n = 6, 3 and 2 for 
February, March and April, respectively). Bars ± 95% CI. 



The gene expression level of HSP70 and P53 showed an ascending 
increase from February to May, with a peak in March/ April. These 
trends mimic the temporal variation in UV recorded between 
February and April for the Gulf of California suggesting that over- 
expression of repair genes is dose- dependent in whales as is known to 
occur in humans and laboratory animals^^'^^. It is reasonable to hypo- 
thesize that the comparatively lower levels of expression recorded in 
May are due to acclimatization to UV exposure, as UV levels pla- 
teaued in April. This phenomenon has been described in humans, 
whose sensitivity to sunburn decreases with increasing frequency 
and duration of solar exposure^°. 



Blue whales appeared able to increase melanocyte abundance in 
response to the increase in UV from February to May, suggesting that 
individuals of this species are capable of modulating their degree of 
pigmentation. Although statistically non-significant, melanin den- 
sity appeared to follow the same increasing trend over months 
suggesting tanning ability (Table 84; Fig. 3C). Tanning has already 
been described for some wild species. For instance, hammerhead 
sharks gradually increase integument melanin in response to an 
increase in direct UV exposure^^ Other fish species^^ as well as fresh- 
water zooplankton^^, have also demonstrated sun-tanning capability. 
Interestingly in this concept, we did not observe analogous changes 
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Figure 4 | Seasonal variation of levels of mitochondrial DNA (mtDNA) 
lesions (A) and melanin density (B) of blue whales. Season day indicates 
the moment the biopsy was taken (day 0 through 20 = February, 21 
through 50 = March, 51 through 80 = April and more than 81 = May and 
June). Only individuals biopsied for the first time in that season were 
included. Each dot represents a different individual. Lines and shaded areas 
represent the gam regressions and their corresponding 95% CI, 
respectively. 

in the pigmentation of the comparatively darker fm whale. It is 
possible that fm whale constitutive skin pigmentation is sufficient 
to counteract the harmful effects of UV. Indeed, fin whales have 
recorded the lowest prevalence of sunburn lesions compared to blue 
and sperm whales^. 

The observed differences in pigmentation plasticity amongst 
species might entail their distinct migratory behaviour. The higher 
prevalence of microscopic lesions^^ and mtDNA lesions at the begin- 
ning of the season might reflect the time needed for UV acclimatiza- 
tion^°, especially as blue whale pigmentation seems to then increase 
gradually. Unexpectedly, towards the end of the season, blue whale 
melanin density appeared to diminish while mtDNA lesions peaked 
once more. Although statistically non-significant, it is tempting to 
speculate that the trend might reflect skin turnover, or cell survival, 
linked to a focal overload of genetic and cellular skin damage within 
an overall field effect of UV-exposure^^'^^. On the whole, the observed 
trends suggest that blue whales are able to modulate pigmentation in 
response to the fluctuating levels of UV through the season, suggest- 
ing marked phenotypic plasticity. 

Despite the opportunity to study these long-lived oceanic preda- 
tors over decades, our findings demonstrate how poorly we under- 
stand the impact of UV exposure and the basic processes that are 
engaged by this environmental insult. The discovery of an apparent 
plastic pigmentation response as well as the use of distinct strategies 
to counteract harmful exposure to UV amongst whale species raise 
questions about the selective pressure that sun exposure has exerted 
on these marine mammals. 

Methods 

Ethics statement. Sampling was approved by the Ethics Committee Regulations of 
the Zoological Society of London and complied with national and international 
regulations on animal welfare. 

Whale sample collection. Cetacean surveys and sample collection (106, 55 and 23 
skin biopsies of blue, fin and sperm whales respectively) were conducted in the Gulf of 
California, Mexico, between January and June 2007-2009 as was described in detail 
previously*^. Resampled individuals were excluded from the analyses, the first sample 
obtained per individual being the one included. Minimum age was estimated by the 
age at the first sighting of a particular individual using the CICIMAR sighting history 
of photo -identified blue whale's data set. 

Skin pigmentation indices using histology analysis. Melanocyte density, used as a 
surrogate measure of constitutive pigmentation*^'^^ and melanin abundance were 
calculated for each individual within a standardized area using digital photographs of 



H&E stained skin sections. Histological and melanocyte density quantification 
methods were performed as described previously*^. Melanin abundance was measured 
using the image processing program Image J (for details see Fig. S2 and Fig. S3) . Tissue 
analyses were conducted on all samples that were unaffected by technical issues such 
as tissue fixation or section size. Details on sample numbers are shown in each figure 
legend or table caption. 

Gene expression assays. Gene expression assays were described in detail previously^''. 
Briefly, total RNA was extracted using the RNeasy® Mini Kit (Qiagen, UK). 
Complementary DNA (cDNA) was obtained by reverse transcription using the 
QuantiTect® Reverse Transcription Kit (Qiagen, UK). All qPCRs were performed in a 
7300 Real-Time PGR System (Applied Biosystems, UK) using Power SYBR Green 
PGR Master Mix (Applied Biosystems, UK). As internal control genes we used the 
genes coding for the ribosomal proteins S18 {RSP18) and ribosomal proteins L4 
{RPL4), which were selected as best intra- and interspecies control genes among four 
genes tested using the packages BestKeeper, geNorm and NormFinder^^. Target gene 
{HSP70, P53, KIN and TYR) primers were designed using sequences obtained in the 
NCBI GenBank database (for sequence details of the primers and corresponding 
targeted gene sections see Table S5). In total, we obtained high quality gene expression 
data for 60 samples including 22 blue whales, 22 fin whales and 16 sperm whales. 

Mitochondrial DNA lesions. To isolate DNA for assessment of mtDNA damage, we 
used phenol chloroform extraction followed by ethanol precipitation. UV-induced 
mitochondrial DNA (mtDNA) lesions were detected and quantified using 
quantitative real-time PGR (qPGR) as described previously^^. The principle of the 
assay resides, under quantitative conditions, in the fact that damaged mtDNA amplify 
at a lower efficiency rate than undamaged mtDNA and thus the threshold crossing 
point (Ct) value obtained is directly proportional to the level of damage^^. 

Statistical analysis. Gene expression levels were analysed using the relative 
quantification method (level of expression of the target gene relative to internal 
control genes) that is based on the ACt method (Cttarget gene - geometric mean Ctcontroi 
genes) ^"^^ As lowcr ACt valucs represent higher levels of expression, it was easier for the 
interpretation of the results to negatively transform these values. Linear models were 
constructed to investigate interspecific, intraspecific and temporal variation in levels 
of pigmentation and mtDNA lesion abundance. For gene expression analyses, linear 
mixed effect models were constructed^'^'^^ Models were built in R and we used a top- 
down strategy to determine which variables explained a significant fraction of the data 
(see supplementary information for more details) Any violation of normality or 
homoscedasticity assumption was corrected by logarithmic transformation of the 
response variable. 
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